Compensation, interstitial defects and ferromagnetism in diluted semiconductors 
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We present a quantitative theory for ferromagnetism in diluted III-V ferromagnetic semi- 
conductors in the presence of the two types of defects commonly supposed to be responsible for 
compensation: As anti-sites and Mn interstitials. In each case we reduce the description to that 
of an effective random Heisenberg model with exchange integrals between active magnetic impu- 
rities provided by ab initio calculation. The effective magnetic Hamiltonian is then solved by a 
semi-analytical method (locally self-consistent RPA), where disorder is treated exactly. Measured 
Curie temperatures are shown to be inconsistent with the hypothesis that As anti-sites provide 
the dominant mechanism for compensation. In contrast, if we assume that Mn interstitials are 
the main source for compensation, we obtain a very good agreement between the calculated Curie 
temperature and the measured values, in both as-grown and annealed samples. 



Diluted Magnetic Semiconductors are materials where 
the interplay of transport and magnetic properties open 
the perspectives of exciting applications. The III-V semi- 
conductors are particularly promising since a low con- 
centration of magnetic dopants can give relatively high 
Curie temperatures for ferromagnetism^ Q. In these 
materials it is found that the Curie temperatures depend 
strongly on methods of preparation and sample history: 
for the same nominal concentration of magnetic ions the 
Curie temperature may vary by large factors. System- 
atic studies show that different annealing treatments dis- 
play a clear correlation between the Curie temperature 
and the conductivity. This indicates that the process of 
magnetic doping is more complex than a straight sub- 
stitution (Mn(Ga)) of (formally) Ga^+ sites by Mu^+ 
atoms, providing a localized magnetic moment and an 
itinerant hole. In fact the original samples are "compen- 
sated" , i.e. the density of holes measured by transport 
is lower than the concentration of magnetic ions due to 
additional donor impurities especially in the samples as 
grown by Molecular Beam Epitaxy (MBE) . The increase 
in Tc after annealing is then interpreted as removal of 
the defects, resulting in an increase in the hole concentra- 
tion which mediates the magnetic exchange. This leaves 
obscure the precise form of the compensating defect, and 
does not provide a quantitative theory relating the Mn^+ 
concentration, the hole density and density of compen- 
sating defects to ferromagnetism. 

There are two probable candidates for compensation: 
both Arsenic antisites Asca (i.e. As atoms on sites of 
the Ga sub-lattice) and Mn interstitials Mn/ have long 
been known to be double donors. The two forms of de- 
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fects differ in an important way: for each AsQa there 
are two holes removed, i.e. only the carrier density is 
changed, while each interstitial, in addition, introduces a 
magnetic moment, changing the number of magnetically 
active ions. Microscopic calculations indicate that the 
Mn/ are preferentially situated on interstitial sites adja- 
cent to occupied Mn(Ga) and that the coupling between 
interstitials and the adjacent moment is essentially given 
by antiferromagnetic superexchange coupling ( J « -320 
K) 3] . In fact there are two inequivalent interstitial posi- 
tions, the Mn atom can be located inside the tetrahedron 
formed by either four Ga T(Ga4) or else four As T(As4). 
We shall make no distinction between the two possible 
positions sites, since in either case the AF exchange with 
Mn(Ga) is strongly antiferromagnetic, and we refer to 
Mn/ as the sum of the two. 

The immediate question is, what is the proportion of 
the two defects: interstitials Mn/ and Asca in the ferro- 
magnetic samples? A conclusion of this paper is that 
the observed Curie temperatures in samples at differ- 
ent stages of annealing, can only be explained assuming 
that interstitial defects dominate cornpensation. Such a 
dominance agrees with Wolos et al. who estimated, 
from the strength of optical transitions, a relatively small 
(fewer than 10% of the total Manganese atoms) num- 
ber of antisites AsGa7 fewer than 10% of the total Man- 
ganese atoms, and, from Electron Paramagnetic Reso- 
nance, a much larger number of other compensating de- 
fects. Similarly Wang et al showed that the saturated 
magnetization at low temperatures was consistent with 
the elimination of two magnetic moments with each im- 
purity. Furthermore polarised neutrons reflectometry 61 
and Auger spectroscopy and resistivity measurements 21 
showed that the annealing process corresponds to redis- 
tribution of Mn sites and the increase of the magneti- 
sation far from the surface. We emphasise that clear 
proof of the role of interstitials is still necessary, as other 
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techniques, by Transmission Electron Micrography 0, or 
by Infrared Absorption and Positron Annihilation Spec- 
troscopy suggested a much higher concentration of 
anti-sites. The new element we are bringing here, is a 
quantitative theory for the Curie temperature, which as 
we shall explain below, is much more accurate than Zener 
mean-field theory. We note that in ref we anticipated 
the fact that the changes in the carrier density due solely 
to anti-sites were insufficient to explain the reduction of 
Tc- 

Recently, by combining first principles calculations and 
an semi-analytical approach, we were able to provide an 
excellent agreement between the calculated Curie tem- 
peratures and those measured in optimally doped semi- 
conductors . In the first step of this method, we derive 
the exchange integrals between magnetic impurities us- 
ing the Local Density Approximation (TB-LMTO) and 
magnetic force theorem |lOl| providing an effective classi- 
cal random Heisenberg Hamiltonian. Note that the cal- 
culated exchange integrals, include the effect of disor- 
der within a Coherent Potential Approximation (CPA) 
for electronic motion. In the second step, we treat the 
random effective Heisenberg model within an approach 
where thermal fluctuations are treated within a local ran- 
dom phase appoximation (RPA), whilst the disorder is 
treated ^exactly\ i.e. the magnetic properties are calcu- 
lated for individual configuration of disorder generated 
by sampling techniques. This theory is an extension of 
ref.|lll| where disorder in the effective Hamiltonian was 
treated by a form of CPA. An attractive feature of an 
"exact" treatment of disorder is that it allows us, for 
exainple, to study the effect of correlations in the disor- 
der [13 , a question of importance in interpreting experi- 
ments on high Curie-temperature samples grown using 
the OMVPE technique (Organo-Metallic Vapor Phase 
Epitaxy) [13 . We attribute the success of our approach 
to (i) the realistic calculations of the exchange integrals 
(ii) to a proper treatment of the thermal fluctuations and 
disorder of the effective Heisenberg Hamiltonian. This 
second point is confirmed by consistency with Monte- 
Carlo simulations 0| . 

We first analyze the dependence of the Curie tem- 
perature Tc with the density of As anti-sites j/^s- 
in our preliminary study of this issueQ, which was re- 
stricted to a single nominal concentration of Mn, we do 
this by introducing in the ab-initio stage of our calcu- 
lation a concentration of anti-sites, which influences the 
calculated exchange couplings. This calculation also al- 
lows us to verify, within LDA that there is indeed the 
compensation expected. In Fig. 1 we plot, for differ- 
ent concentration XMn of Mn, the variation of the calcu- 
lated Tc in (Gai_a;^j„_y^^Mnj:„_^Asy^JAs as a function 
of 7 = Uh/xMn- Since each As anti-site is a double donor, 
the carrier density is Uh — xun — '^DKs- We observe that 
above a critical value 7c(a;Mn) , Tc is weakly sensitive 
to As-antisites, this is particularly clear for xmh = 0.03 
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and 0.05. For 7% Mn, for example, we also observe that 
for 7 < 0.50 ferromagnetism is unstable. The reason for 
this is that as the density decreases, the nearest neigh- 
bour exchange integral becomes increasingly dominated 
by the (antiferromagnetic) superexchange contribution, 
leading to frustration. Note also that this behaviour with 

hole density is incompatible with Zener mean-field theory 
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which predicts that Tc oc . Let us now discuss the 
compatibility of these results with the assumption that 
As anti-sites dominate the mechanism of compensation. 
We do this by comparing samples with approximately 
fixed total density of Mn impurities but exhibiting large 
variation in their Curie temperatures. 

We plot in Fig. 2, the measured as a function 

of 7 = nh/xMn- In contrast to calculated values, Tq^^ 
is more sensitive to the carrier density, varying linearly 
with 7. The other important difference with the curves 
in Fig.l is that ferromagnetism is still observed for rather 
small values of 7. Thus, if we assume that As anti-sites 
dominate compensation, theory and experiment would 
disagree. As our approach was successful for uncompen- 
sated samples, and consistent with Monte-Carlo [l^ we 
conclude that As anti-sites do not dominate compensa- 
tion. 

As already mentioned, the saturated magnetization[^ 
at low temperature indicates that the compensating de- 
fects affect both the density of carriers and the density 
of magnetically active Mn impurities. Let us now take 
the alternate limit in which compensation is taken to 
be entirely due to the presence of Mn interstitial defects 
Mn/. We denote by xun , ^CMnoa ^.nd XMn(I) respectively, 
the total density of Mn, the density of Mn on Ga sub- 
lattice and on interstitial location respectively. Recent 
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FIG. 2: Experimental values of Tc (,2]) as a function of the 
measured hole density per magnetic impurity 7 — for 
Gai-2,Mn2,As at nominal concentration 6.7%. 



FIG. 4: Tc for Mna;Gai_3; As. The experimental data (di- 
amonds) are plotted assuming for each sample a;Mn = x^s 
and 7 = 7cff (see text). The values of 7eff are shown in the 
figure. The small circles , the square and triangles are cal- 
culated Curie temperature. The small circles (dashed line) 
correspond to uncompensated samples 7 = 1 and the square 
and triangle correspond to xiAn ~ 0.035 and 0.05 respectively. 
For these cases, the value of 7 is shown in the figure. 



t ^ if ^ 



Effective probiem with: 



• Mn(Ga) 

• Mn(l) 
O hole 



X eff=xtot-2x(i) 
n h=xtot-3x(l) 
xtot=x(Ga)+x(i) 



FIG. 3: Left side: the up (resp. down) arrows indicate spin of 
Mn(Ga) (resp. Mn/). The small circles are itinerant carrier 
(holes). Mn/ are double donor and strongly coupled antifer- 
romagnetically to Mnca- Right side: effective model with XeS 
MuGa impurities and Uh holes. 



first principle calculations [3 and channeling Ruther- 
ford backscattering experiments indicates that Mn 
interstitials are preferably attracted by MuQa and tend 
to form pairs of spins with strongly antiferromagnetic 
couplings. Thus, we suppose that Mn/ are not com- 
pletely random, but are only in positions with a Mnca 
as a nearest neighbour (see fig. 3). In writing an ef- 
fective Hamiltonian, we will eliminate the strongly anti- 
ferromagnetically coupled pairs of the Mn/ and adjacent 
MuQa- They can be assumed, within high precision, to 
form bound singlet pairs whose effect on the magneti- 



cally active ions is small [l^. The remaining "active" 
Mn of density XeS = xun — 2a;Mn(I) which are not di- 
rectly coupled to a Mn/, and with the measured carrier 
density n^, interact via an effective Heisenberg model 
with couplings determined by the measured carrier den- 
sity rih- We make the same calculation as before but with 
the effective concentration XcB and compensation 7e//- 
Since each Mn/ is a double donor, the total density of 
carriers is Uh = a^MnCGa) - 2a;Mn(I) = xun - 3a;Mn(I)- 
Thus from each measured we can deduce the density 
of Mn/ and of 'unpaired' Mn which are, respectively, 
2:Mn(I) = l{xun-nh) and XcB = i(a;Mn + 2n,i). We also 
define the effective 7 parameter as 7eff = = — — . 

In Fig. 4 we show both the experimental data and cal- 
culated Curie temperature for various 7 as a function of 
a^Mn = XcS- The variable XgS are calculated from the 
values given in and the corresponding values of 7efi is 
noted on the figure for each sample. We take the xun 
to be the nominal concentration of each sample; thus 
we neglect effects of surface inhomogeneities^ 0, Q in 
calculating the bulk ordering temperature. First, we ob- 
serve that the well annealed samples (of highest Tc) are 
in excellent agreement with the calculated values for un- 
compensated samples, 7 = 1 curve ('optimal curve'). We 
remark that this optimal curve, which depends on ex- 
change integrals recalculated for each concentration, can 
nevertheless be well parameterized by the simple form, up 
to XMn = 0.10, Tc « A{xMn " Xc)'^/'^ whcrc Xc = 0.0088, 
A = 649 K for Ga(Mn)As. The samples corresponding 
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to intermediate Tc are still in reasonable agreement with 
the optimal curve for uncompensated samples at the ef- 
fective concentration of magnetically active ions. The 
deviation from this optimal line is small but increasingly 
visible for as-grown samples. In order to refine our calcu- 
lations we have taken into account that for these samples 
7cff is substantially less than 1. We therefore performed 
additional calculations for fixed xmh — 0.035 and 0.05 
and various hole densities. As in Fig. 1, we vary ^^s in 
this effective Hamiltonian with the addition of anti-sites, 
here used as a purely calculational device to change the 
carrier density, whilst keeping the calculation fully self- 
consistent. We now observe that the agreement with the 
experimental measurements is very good for all the mea- 
sured samples (as-grown and annealed). For example the 
as-grown sample which corresponds to Xeff ~ 0.035 and 
7eff ~ 0.52 agrees very well with the calculated value 
(square symbol) for the same parameters. Note that us- 
ing the above relations we find that the density of Mn/ 
in as-grown samples is a:Mn(I) ~ 0.016 which corresponds 
to approximately 25% of the total Mn density. This is in 
very good agreement with the value estimated in ref . [l7| 
(see Fig.4 of ref.firj). 

In conclusion, we have shown that experimental mea- 
surements in samples with fixed nominal magnetic impu- 
rity concentration could not be explained assuming As 
anti-sites as the dominant mechanism for compensation. 
On the other hand, with the assumption that Mn in- 
terstitials dominate, we obtained an excellent quantita- 
tive agreement with the measured Tc in both as-grown 
and annealed samples. It may be possible, in varying 
sample preparation, to increase the number of anti-sites, 
this will have a weak effect on the Tc, provided the 
7e// remains above the region of instability, as seen in 
Figure 1. For Ga(Mn)As samples, we can write an ex- 
plicit first approximation (7cft — 1) using the empirical 
form for the optimal curve, by replacing XMn by x^s'- 

Tc W ^( ^M„+2n;, _ 0.0088)^/2 ^^^^^ ^ ^ g4g rp^ 

take into account the smaller corrections due to the value 
of 7cff we do not have an explicit analytical form, but 
numerical corrections can be predicted as in Figure 1. 
Hence our combined ab-initio/local-RPA approach is a 
very powerful tool to study ferromagnetism in diluted 
ferromagnetic systems even in the presence of compen- 
sating defects. The same approach can be applied to 
macroscopic inhomogeneities, for example surface effects, 
which may be necessary to understand thin films and de- 
vices. 
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